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We present an approach to solve a Bethe-Salpeter (BS) equation exactly without any approxima- 
tion if the kernel of the BS equation exactly is instantaneous, and take positronium as an example to 
illustrate the general features of the solutions. As a middle stage, a set of coupled and self-consistent 
integration equations for a few scalar functions can be equivalently derived from the BS equation al- 
ways, which are solvable accurately. For positronium, precise corrections to those of the Schrodinger 
equation in order v (relative velocity) in eigenfunctions, in order v in eigenvalues, and the possible 
mixing, such as that between S (P) and D (F) components in J^ = 1 (J^ = 2"*"+) states as 
well, are determined quantitatively. Moreover, we also point out that there is a problematic step 
in the classical derivation which was proposed first by E.E. Salpeter. Finally, we emphasize that 
for the effective theories (such as NRQED and NRQCD etc) we should pay great attention on the 
corrections indicated by the exact solutions. 

PACS numbers: ll.lO.St, 36.10.Dr, 12.20.Ds 

Keywords: instantaneous BS equation, exact solutions, positronium, relativistic corrections. 



Bethe-Salpeter (BS) equation |l| is a very good tool 
to treat various bound state systems. For a fermion- 
antifermion binding system, the BS equation is written 
as follows, 



'nii)xp{q){h + m2) = i 






ViP,q,k)xp{k) , 

(1) 

where Xpil) is the BS wave function, P is the total mo- 
mentum, q is relative momentum, and V{P, k, q) is the 
so-called BS kernel between the electron and positron 
in the bound state, pi,P2 are the momenta of the con- 
stituent electron 1 and positron 2, respectively. The to- 
tal momentum P and the relative momentum q are re- 
lated to the momenta of the two quarks as follows: pi = 
aiP + q, Qi = "1 , P2^a2P-q, a2 = "" . 

^ ^' ^ rni+m2 ' -^ ^j .^ m\+ra2 

If the kernel of the four-dimensional BS equation, 
V{P,k,q), in P = frame (C.M.S) of the concerned 
bound state, has the behavior: V{P, q, k)\p^g — V{q, k) , 
the BS equation is called as 'instantaneous one', and may 
be derived to a Schrodinger equation accordingly, that 
was firstly realized by E.E. Salpeter [21 . 

In Coulomb gauge, the transverse photon exchange 
is considered as higher order, so at the lowest or- 
der the BS equation for positronium has the kernel 
V{P, q, k)\postromum = 7°K7" = "7° 7^7° Only, that 

y > 

is instantaneous, so we take positronium as an example 
to pursue exact solutions of the BS equation. It certainly 
is interesting that the relativistic corrections, including 
the possible mixing, such as that between S and D com- 
ponents in J^'-^ = 1 state and that between P and F 



components in J^'~^ = 2++ state etc, will be fully fixed 
by the instantaneous BS equation. Here we report that 
as the case of positronium, an instantaneous BS equation 
can be really solved without any approximation, and out- 
line the approach briefly. We derive the instantaneous 
BS equation (the BS wave function in 4 spin structure) 
into a set of coupled and self-consistent integration equa- 
tions for its components (scalar functions) without any 
approximation, and solve them numerically. The accu- 
racy of the solutions can be contral at all. The results are 
discussed finally. While the details are put in Refs.l6l,'3|. 

The approach which we show here mainly is to follow 
E.E. Salpeter derivation but without any approximation. 

As done in Ref.^ and in book ^, the instantaneous 
BS wave function (pp{q) is introduce: 



</5p(q) 



"2^ 



xp{q°,q), 



then the BS equation Eq^is re- written 



(2) 



Xp{q",q)^Sf\p'iMq)Sfn.~p'2) 



(2) 



(3) 



where S } (pi 



andS'l^^(- 



r yjjij o,i±ii ^ ,■ \ P2) are the propagators of the 
fermion and anti-fermion respectively and the integrated 
'BS-nut' ri[q) = / j^^^^^ k)tpp{k) . For general applica- 
tions, we keep toi 7^ 7712 at this moment, although final 
application in this paper is to positronium mi = 7712. 
The propagators can be decomposed as: 



A^ 



iJS)'>{Jp^^) 



t^\ _ 



Ki9) 



K («1 



JqO+UiM-uJi+ie ' .Jq°+aiM+LJi-ie ' 



(4) 
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with uji = y/mf + (f , ^f{q)^^^ l°^i ± J {mi + 7 • 
q) , where J ~ 1 for the quark(i = 1) and J = — 1 for 



the anti-quark(z = 2). It is easy to check 

Af{q)j^Af{q)=0, Afiq)j°Af{q)=Afiq), (5) 

Namely A* are 'energy' projection operators and com- 
plete. With them for below discussions let us define 

ipp^iq) as: 



Oa±/ 



^^^iq)^Af{q)^\piq)^^Afiq) 



(6) 



Because of the completeness of the projectors A* , for the 
BS wave function ipp{q), we have: (pp{q) = (fj, (q) + 
ipp^{q) + (pp'^{q) + ipp^{q) . As done by E.E. Salpeter, 
to derive the instantaneous BS equation to corresponding 
Schrodinger equation, by carrying out a contour integra- 
tion of the time-component q'^ on both sides of Eq.©, 
we obtain 

,^ (^ _ A+(.r)»;p(g-)A+(gl _ A^iq)npiir)^^{q) ^7^ 

fFVl) {M-UJ1-UJ2) {M+UJ1+UJ2) ' ^ ' 

and applying the complete set of the projection oper- 
ators Afp{q) to Eq.O further, we obtain the coupled 
equations: 

(M -uJi- u:2)^+p+{q) = A+ {q),jp{q)A+ (q) , (8) 



{M + uJi+uj2)^p-{q) - -A^{q)rjp{q)A^{q), (9) 



<y5p io) ^^p'^io) = 0. 



(10) 



Note that because of the completeness of the projector 
set, Eqs.(|Hl El [Tn|) in whole are exactly equivalent to 
Eq.0 and Eq.Q. 

While E.E. Salpeter |3| and the authors of 3] would 
like to connect the above coupled equations to the Breit 
equation |j,|a|, so they 'combined' Eas. H8l9llO() into one 
operator equation (in C.M.S. of the bound state): 

[M ^ Hi{q) - H2iq)Mq) 
- A+iq)j"7^{q)^^A+iq) - A" (9)7^(9)7° A-(g) , (11) 



with the definitions Hi{q) = miP + q- a , i?2(<?) = ni2l3 — 
q ■ d , P — 7° , a = Pj. Namely they considered the 
equation Ea. (|ll|l with the definition i?i.2 as an operator 
representation of the coupled equations Eqs. (|S1 El I10|l . In 
fact, it is no correct. Ea. ljlll) is not fully equivalent 
to the coupled-equations Eas.(pnr^ll()|l. 

Now let us show the un-equi valence as follows. 

When applying the project operator A'^{q)j'^ (g) 
7°Aj((j) to Ea. ((TT|l we obtain Eq.®, when applying the 
project operator A'^ {q)j^(^j^A2 (q) to Ea. ((TT|) we obtain 
the Eq.Q, whereas when applying A^(q)^^ (g) 7" A J ((7) 
to Ea. Hmi . then we obtain: 



[A/- 



■uji{q} - 
-uji{q) 



-^2{q)]v^' 
-uj2{q)](p^' 






0, 
^0, 



with iLJi2 — \ nrii 2 + 9^- When mi = 7712 i.e. uii = 1^2, 
we may obtain Eq. pU|) from Eq. (|12|) unless M = in the 
extremely relativistic cases, thus in the cases mi = 7712, 
equations Ea. H12|l are equivalent to Ea. pOfl . Whereas 
when TTii 7^ 7772 i-e. wi 9^ W2, the equations Ea. l|12|l 
have the 'trivial solutions' just as Eq.lQUj), but also have 
'non-trivial solutions'. Practically when solving equa- 
tions, such as Ea. (|ll|l . one would not check the pro- 
jections A+(g)7" (g) 70Aj(g), Ai{q)-i° (g) 7°A^(g) and 
Aj''((j)7° (g 7° A J ((7) on the solutions precisely, and gen- 
erally non-trivial solutions are chosen. Especially, when 
further approximations are made, certain misleading may 
occur i.e. the obtained solutions may not satisfy the 
equations Ea. H10|l . even in the cases 7771 — m.2 6J. There- 
fore, we conclude that Ea. H12(l are not fully equivalent to 
Eq.Jini). 

We will present the different consequences due to 
Ea. (|10() and Eq.l|12|) more precisely and solve them re- 
spectively in Ref. J] . 

In the literature the authors of Refs.0, Q etc made 
an additional assumption on the spin structure of the 
BS wave function e.g. the spin structure of the wave 
functions for 5- wave positronium has the formulation (in 
C.M.S. P = 0): 

'^0- + (i5o) ( 9 ) = [70 + 1] 75/(9^) (13) 

for JP^ = {)-+ (i^o); and 

<Pi-(3Si).a(9) = [70 + AAx.f{q) (14) 



^a7+-£^^. 



for J^^ = 1— (^S'l). Here Ax = e+r 
With the spin structure finally the bound state problem 
turned to solve a Schrodinger equation accordingly by 
taking the so-called positive-energy equation Eq. Q only 
(the rests are ignored). 

To solve the instantaneous BS equation Eq.© for 
positronium without any approximation, even in the case 
7771 = 7772 we 'retum' to solve the Eas.lpi^ lTUI) . In fact, 
the approach here is also applicable in the cases 777,1 7^ '772. 

In P — frame (C.M.S.), the most general formulation 
of the BS wave function for the bound state J^'~^ = ^ 
(^5*0) may be written as the follows: 

</'iSo(9) =7% VdgD +7^2(191) 

V 3 2777 




37— i^^4(M),(15) 



where E 



^1-17"'" + ^117 - ^107' 



1' 



1 +il 
V2 



7 = ^ ^"^ and Yj^ = Yhn{Oq,(f)q) are spherical har- 
monics. To apply the equation Ea. llU|) to Ea. (|15|) . we 
obtain the 'constraints' for the components of the wave 
function precisely: 



(12) 



1 



Mm^-i,Mm, Mm = o 



(16) 



TABLE I: The expansion coefficients C^''^, for the J^'^ 
states. 



NL 


WF 








^iAS 


(I'So) 


/r' 


-0.707104 


8.934D-6 


4.225D-6 


2.629D-6 


/^^' 


-0.707109 


6.962D-6 


3.206D-6 


1.979D-6 


(2^50) 


/}^' 


5.031D-6 


0.707106 


-7.246D-6 


-3.571D-6 


I2 


7.004D-6 


0.707107 


-6.567D-6 


-3.169D-6 


(3^50) 


/f' 


2.188D-6 


5.898D-6 


0.707107 


-6.768D-6 


j-(3) 
I2 


3.207D-6 


6.577D-6 


0.707107 


-6.425D-6 



for J^'-^ — +(^S'o) states, in fact, can be simplified to 
write as 

^^M - (/ + 1 + y^i?/^|/7^ (20) 

I V "J ^T^ I 



with /i(-Z) ~ /2(9) = / = Rn(l=oM)- 

Let us now solve the equations for the 1 states. 

Similar to '" states, we start with the most general 
BS wave function, and have the constraints Ea. 1)101) ap- 
plied, then the BS wave function becomes the form: 



(^3, (^4 in Ea. (fT5)l are replaced with Ea. (fT?)|l . then the rel- 
ativistic wave function of the state '" contain two inde- 
pendent components (^1,(^2 only. With straightforward 
calculations, from Eqs.|(Slini fo'' positronium, finally we 
obtain the two coupled radius equations: 

(M - 2m) /i -I- 2m/i - 2m/2 = 



0.3 zn 



/id|fc|Qo/2, 



(M - 2to)/2 - 2^/1 + 2m/2 = 



(17) 



where Q„ = Q»( ^m g ) ('^ = 0' l^ ' ' ') ^'^'^ the n the 
Legendre functions of the second kind. Here to shorten 
the notations we have introduced fi^2'- '^i,2{\<f\) = /1.2 
and Lo — wi^2, and the angular integrations have been 
carried out precisely. Now without any approximation 
a self-consistent problem, two coupled equations Eq. 1)17)1 
for two independent scalar functions /i,2, is reached to. 

Then we further solve the coupled equations Ea. lfT7|) 
numerically by expanding the components /i 2 and the 
equations in terms of the bases of the exact solutions of 
the Schrodinger equation [9j, and then diagonalizing the 
matrix equations. 

The final results for the eigenvalues E'^^^ , j = 1, 2, • • •: 



^(1) =-6.8026952534; E^^^ 
^(3) = -0.75586715480; E^'^'i = 



-1.70069524809; 
-0.42517586843; (18) 



and the corresponding eigenfunctions 



frm)^Y.^ij.i-Rnm), 



(19) 



<^i--,a(9)= 2yfCA + ^A /l 



2x/¥^A7' + ^A7" /2 



2.,/fCA + Sj/3 



+ 



Dxl^ - i?A7° /4 , 



(21) 



where Ax is the same as in Eg . 1)1411 and 

Bx = V6e+-/+Y2-2 - V3(^4l^ + e^7+)>"2-i 

+ (6+7- +6^7+ + 26^7^)^20 

~V3(e-j^ + efj-)Y2i + V6e;^7">^22 , 



Cx = 



\q\ 



4Y1-1 + ^xYn - 4Yio 



Dx^^4 



,-^A^ .A_ 



l-)Yn + (etl+-4l'^)Y,^i 



+ (47"-eA7+)no 



Here A = is introduced for the 0*'' component to avoid 
confusions. From Ea. l)21)l . we may see clearly that 73^4 
correspond to _D-wave components while /i_2 to 5- wave 
ones. The wave function now has four independent com- 
ponents /i(i = 1, . . . ,4). 

From Egs. 1)819)) and to carry out the angular integra- 
tion in the equations, finally we obtain the coupled radius 
equations: 

(M - 2m) f, =. |^¥/4 - 2mA - ^^^^A 



-f/f^N 



^mn,;,-(|mQ, 



Qtt muj 



}Qo)f: 



nl 



where (j) denotes those correspond to the jth eigenvalue 
Eq.ldHl, Rnii\q\) (nl = IS, 25, 2P, 3S, 3P, • • •) are the 
Schrodinger radius solutions in momentum representa- 
tion for the positronium and their precise expression can 
be found in Ref.Q- For the low-laying solutions the 
coefficients Cp2; (nl < 45) are put in TABLE I and 

Cl^rd < W-^ {nl > 45). Based on the values in TA- 
BLE I, up to the order v accuracy the eigenfunctions 



_ 2^5 ^ f^ _ 22vf+u^f^ _ 2mf^ 



(A^ -2m)/2 = 01^/3- 



'7/f¥l 



5 |g1|fc|, 



1/3 - {\ 



MIM, 

3 inuj 



h + ^Qo /i 



(M - 2m)/3 = -2m/3 -f f^^^A - W€ h 



■ / %d\k\ [ (|1|q, + ^Q,) /, - I^MMq,/, 



3V5 

^J 

3V5 mw 



(Af - 2m)/4 = i2l»!±4^/3 _ 2m/4 - ^M!/^ 



■f/§^l^l[( 



2k1|fc|( 

3 muj 



3V5 "I 
1 ) fs 



2y^\q\\k\ 
3V5 "li^ 



Qi/i 



(22) 



Four coupled equations for four independent components 
/i,2,3,4 is the request of the S-D wave mixing and a self- 
consistent problem is reached. 

We solve the Eas.()22ll numerically with the same 
method as adopted in the case for ^, and obtain the 
results: 

£;(i'(l3S'i) = -6.8027275 ; E'^^\2^Si) = -1.7007008 ; 
£:(3)(335'^) ^ -0.7558690; E'^^^S^Di) = -0.7558729; 

(23) 

for eigenvalues, while the table of the coefficients C^ „; 
for the eigenfunctions is too big to present here. Instead, 
we present them here only up to the order v accuracy as 
in the case for ^. 

For J^^ = 1 (^Si) states we can simplify: /i(?) ~ 

-/2(g) = / = i?„(,=o)(|g1) and h(q) ^ fiiQ) ^ 0, i.e. 



f^SuxC^) 



V 3 m '^ 



A: 



2\/i^^A7' + ^A7^ 



/. 



(24) 



For J^'^ = 1 (■^-Di) states we can simpHfy: fi{Q] 
/2(g) ~ and /3(9) ~ /4(9) = / = i?„(z=2)(|g1), i.e. 



^^Di.xil) 



-2 



3 rn 



Cx + B, 



^MZ5^^5_5^^0 



(25) 
can be 



The accurate coefficients C^ J, for J^'-^ = 
iound in our achieved paper |3]. 

We may see that the exact solutions obtained here 
for the wave function contain the order v = — correc- 
tions to the approximate ones obtained by E.E. Salpeter 



approaches (described as above), namely, the terms in 
Eqs.l^OI [^ |23J) explicitly proportional to v. The cor- 
rections will cause the normalization for the wave func- 
tion a change in order w^ (the detail in Q). The most 
interesting fact is the exact solutions for 1 states 
present the possible S ^ D wave mixing. To compare 
the eigenvalues Eq.l^Hl) foi' 0^+ and Eq.ll^HIl for 1 
with those of the Schrodinger solution (the degenerate 
ones without any splitting), we may see the 'hyperfine' 
splitting for spin singlet and triplet states which is in 
order of a^ {v ~ a), and there is a remarkable split- 
ting between npS*!] and n[^Di]{n > 3) for 1 states 
^(3) (335^) > Ei'='' > E^'^^S^Di) and A^ ~ a^ due to 
the fact that the BS equation Eq.Q is relativistic. 

Note that the approach we present here may be ap- 
plicable for all the instantaneous BS equations and can 
obtain the exact solutions, while the features of the exact 
solutions for positronium described above may be differ- 
ent for an instantaneous BS equation with a different 
kernel from that of Coulomb interaction. 

Based on a specific problem of positronium, we have 
learnt that the wave functions contain quite great rel- 
ativistic corrections thus for effective theories, such as 
NRQED 113 (and NRQCD [ij), the relativistic correc- 
tions from the wave functions should be considered care- 
fully if the relevant calculations declare the relativistic 
effects are taken into account. For NRQED (NRQCD) 
etc, the wave function effects may be involved either at 
matching the underlined theory to the effective one or 
into the non-pertubative matrix elements. Moreover the 
problems, such as gauge invariance for the matrix ele- 
ments and El (electrical dipole) radiative transitions etc, 
are sensitive to the effects, the exact solutions are crucial 
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